On the Red-Green-Blue Model 
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We experimentally study the red-green-blue model, which is a sytem of loops obtained by su- 
perimposing three dimer coverings on offset hexagonal lattices. We find that when the boundary 
conditions are "flat", the red- green- blue loops are closely related to SLE4 and double-dimer loops, 
which are the loops formed by superimposing two dimer coverings of the cartesian lattice. But we 
also find that the red-green-blue loops are more tightly nested than the double-dimer loops. We 
also investigate the 2D minimum spanning tree, and find that it is not conformally invariant. 

PACS numbers: 05.50.+q, 64.60.Fr, 64.60.Ak 
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INTRODUCTION 

We investigate the red-green-bluc (RGB) model, which 
was introduced by Benjamini and Schramm. An RGB 
configuration is a system of the loops on a region of the 
triangular lattice, which is obtained by superimposing 
three perfect matchings (or dimer coverings) on offset 
hexagonal lattices as shown in Figure [[} The sites of the 
triangular lattice may be three-colored so that no two 
adjacent sites have the same color. If we delete the sites 
of a given color (say blue), then the sites of the other two 
colors (red and green) form a hexagonal lattice on which 
we can construct a random perfect matching (the blue 
perfect matching). When we superimpose the red, green, 
and blue perfect matchings, each vertex is matched with 
one neighboring vertex of each of the other two colors. 
Since each vertex has degree two, an RGB configuration 
consists of closed loops, and each loop has an orientation 
if we follow the edges in the order red to green to blue. 

It is worth remarking that the boundary conditions of 
dimer systems can have a profound impact on the be- 
havior of dimers even far within the interior of a region 
Jj], ||, U . There are height functions associated with dimer 
configurations on the cartesian lattice || and hexag- 
onal lattice H (see also @, §])■ If there is an imbalance 
between the different colors of vertices along the bound- 
ary, then the height along the boundary will be "tilted" , 
and this affects the dimers throughout the region. Conse- 
quently, the behavior of the RGB model on regions with 
tilted boundary conditions could be different from the be- 
havior of the RGB model on the regions that we consider 
here, where the three color classes along the boundary 
are balanced ( "flat" boundary conditions) . 

In earlier work, Kenyon and the author found ex- 
perimentally that the fractal dimension of these loops is 
3/2. Here we report on additional experiments, where we 
find that the winding angle variance at a typical point 
on a loop is I x log!? where D is the diameter of the 
loop, and that the system of loops appears to be confor- 
mally invariant. These properties suggest that the RGB 
loops are closely related to stochastic Loewner evolution 
Q with parameter k — 4 (SLE 4 ), and that the RGB 
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FIG. f: The red-green-blue (RGB) model on a triangular 
region of side length 12 (L = 4). The three color classes 
of vertices are balanced along the boundary, so the boundary 
conditions are "flat". 



loops belong to the same universality class as the con- 
tours of 2D Fortuin-Kasteleyn jflj] clusters at criticality 
when q — 4, and the loops formed in the double-dimer (or 
"double-domino" ) model p2[ , which in turn are thought 
to correspond to the "contours" of a Gaussian free field 
p"2| . However, the system of RGB loops does not have 
the same limiting behavior as the system of double-dimer 
loops, because we also find that the RGB loops are more 
tightly nested than the double-dimer loops. 
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FIG. 2: The RGB model on the triangular domain (L = 20) is conformally mapped to a disk. The inverse conformal map, 
from the disk to the triangle, is the hypergeometric function shown in the diagram. If the RGB model is conformally invariant, 
then the loops in disk will be rotationally invariant. 



GENERATING RGB CONFIGURATIONS 

To generate an RGB configuration of a region, we need 
to generate dimer coverings of three regions of the hexag- 
onal lattice. There are many ways to generate dimer cov- 
erings of the the hexagonal lattice, but the fastest of these 
is based on a generalization of Temperley's bijection ]Tq| 
between spanning trees and dimers. To generate the per- 
fect matchings on the hexagonal lattice, the correspond- 
ing spanning trees are on a directed triangular lattice (see 
Jl4j for details), and these spanning trees may be quickly 
generated using an algorithm based on loop-erased ran- 
dom walk |ll| . This generalized- Temperley bijection only 
works for regions of the hexagonal lattice that have cer- 
tain special (very flat) boundary conditions, so we can 
only expect to use it to generate RGB configurations of 
certain nice regions. One region where we can use span- 
ning trees to rapidly generate RGB configurations is the 
equilateral triangle with side length 3L (Figures |l| & |^) . 



WINDINESS OF RGB LOOPS 

We recall the definition of the windiness of a loop used 
in (TgJ . Consider an ant which travels along the loop; 
after the ant has just traversed a given edge in the loop, 
before traversing the next edge it will either turn left 
120°, turn right 120°, or not turn at all. If we keep track 
of the total turning (measured in radians) when the ant 
travels from point A to point B on the loop, then this 
is (approximately) the winding angle between points A 



and B. When the ant travels all the way around the loop, 
it has turned ±360°, so to make the winding angle be- 
tween points A and B independent of the number of times 
that the ant travels around the loop and the direction of 
travel, we adjust the total turning (measured in radians) 
by 2ir x [# steps between A and B] 4- [length of loop] . To 
define the winding angle at a given point X relative to 
the global average direction, we pick an arbitrary point 
A, compute the winding angle from A to X, and sub- 
tract a global constant so that the average winding angle 
at points on the curve is 0. 

When we measure the variance in the winding angle at 
random points along the longest loop in the RGB config- 
uration in a region of order L, we find that the variance 
grows like 1 x log L — so in the notation of Jll| , «2 = 1. 
This winding angle variance coefficient of 1 also shows 
up in in the contours of FK clusters at criticality when 
5 = 4, and other related models such as the double-dimcr 
model, and SLE 4 (see @ fj||). 

CONFORMAL INVARIANCE OF RGB LOOPS 

Since there is only one region, the equilateral triangle, 
for which we can rapidly generate RGB configurations, 
this makes the testing of conformal invariance somewhat 
interesting. The test that we use is similar in spirit to the 
tests used by Schramm to test the conformal invariance 
of the uniform spanning tree. We conformally mapped 
the RGB model on the triangular domain to a circular 
domain, as shown in Figure ^. If the RGB model with flat 
boundary conditions were conformally invariant, then it 



must be that after we map a region to the disk, the result- 
ing system of loops would be rotationally invariant. But 
if conformal invariance failed to hold, then there would 
be no particular reason to believe that the loops mapped 
to the disk would be rotationally invariant. After all, re- 
ferring to Figure [^, the points in the disk to which the 
corners of the triangle are mapped certainly look differ- 
ent than other points in the disk, so a priori we would 
expect the image in the disk to be anisotropic if the RGB 
model were not conformally invariant. As we shall see, 
the minimum spanning tree model fails this test, so this 
test is a nontrivial test of conformal invariance. 

To test the rotational invariance of the image of the 
RGB model in the disk, we singled out the outermost 
loop surrounding the center of the circular domain, and 
collected statistics on its furthest extents in the ±x- and 
±y-directions. If the loops in the RGB model are confor- 
mally invariant, then these four random variables would 
be equidistributed. But otherwise, there would be no 
particular reason to believe that any of these random 
variables (other than the first two) would have the same 
distributions. As shown in Figure |3[ the cumulative dis- 
tribution functions for these four random variables ap- 
pear to coincide, so we conclude that the RGB model 
appears to be conformally invariant. 
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FIG. 4: The minimum spanning tree (MST) on the 32 x 32 
square grid, mapped conformally to the disk so that three of 
the corners of the square are mapped to the cube roots of 
unity. The paths connecting these corners is highlighted; the 
triple point T is the point contained in all three paths. 
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FIG. 3: In the RGB model (L = 1024), the outermost loop 
surrounding the center point is selected, and conformally 
mapped to the disk. The image of this loop appears to be 
rotationally invariant, in particular, the empirical cumulative 
distribution functions for furthest extents of the loop in the x- 
direction, {/-direction, — x-direction, and — {/-direction appear 
to coincide. 



CONFORMAL JVOAT-INVARIANCE OF 
MINIMUM SPANNING TREES 



To evaluate the efficacy of our conformal invariance 
test, we applied it to two additional models: the mini- 
mum spanning tree and uniform spanning tree models. 
The minimum spanning tree (MST) is formed by as- 
signing uniformly random edge weights to the edges of 
the cartesian lattice, and picking the spanning tree (con- 
nected acyclic subset of edges) which minimizes the total 



weight. The uniform spanning tree (UST) is simply a 
spanning tree chosen uniformly at random from all span- 
ning trees. 

Figure | shows the MST of a square grid after it is 
conformally mapped to the unit disk, with three of the 
corners mapped to the three cube roots of unity. To test 
the rotational invariance of the MST after it is mapped 
to the disk, we looked at the paths connecting the three 
points at the cube roots of unity (highlighted in Figure ||), 
and focused on the "triple-point" T contained in all three 
paths. If the image of MST in the disk were isotropic, 
then T, e 27ri/3 T, and e 4jrl/3 T would be equidistributed. 
However, as Figure [| illustrates, these variables are not 
equidistributed, so we conclude that the MST is not con- 
formally invariant. The conformal non- invariance of the 
MST is surprising, given the close relationship between 
the MST and invasion percolation jl9| , the close relation- 
ship between invasion percolation and percolation, and 
the conformal invariance of percolation 120, 
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In contrast to the MST, the UST passes this test, as 
shown in Figure The triple point T connecting three 
boundary points of the UST is already known to be con- 
formally invariant p3[ , and indeed the entire UST process 
is now known to be conformally invariant p4j . 

Thus we learn not only that the minimum spanning 
tree is not conformally invariant, but that this test is a 
nontrivial test of conformal invariance. 
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FIG. 5: The distribution of the triple point T for the MST (on 
the 1024 x 1024 grid) is not symmetric under 120° rotations, 
as it would be if the MST were conformally invariant, (u) = 
exp(27n/3).) 
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CONCLUSIONS 

Our experiments indicate that the loops of the RGB 
model (with flat boundary conditions) are conformally 
invariant and have windiness constant 1. Earlier exper- 
iments fl have indicated that the fractal dimension is 
3/2. These properties suggest that RGB loops belong 
to the same universality class as double-dimer loops, the 
fully-packed- loop model with n = 2, and the contours of 
critical FK clusters with q = 4, and are closely related 
to SLE4. However, the system of RGB loops (not just 
individual loops) differs from the system of double-dimer 
loops, because the loops are nested within one another 
more tightly. 
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FIG. 6: In contrast, the distribution of the triple point T 
for the UST (on the 1024 x 1024 grid) passes this conformal 
invariance test. Kenyon has proven that the triple point 
for UST is in fact conformally invariant (see also |24| ) . 



NESTING OF RGB LOOPS 

For a scale-invariant loop model on a region of side 
length L, we would expect the number of loops surround- 
ing a point to scale as const, x logL. The value of this 
constant is a measure of how deeply nested the loops are. 
For the double-dimer model, Kenyon |25| proved that this 
nesting constant is 1/7T 2 . We measured the nesting con- 
stant of the RGB loops, and found that it is 20%-25% 
larger than the double-dimer nesting constant, but we do 
not have a guess for its exact value. Figure shows that 
the outermost red-grcen-blue loops are in a sense larger 
than the outermost double-dimer loops, which is con- 
sistent with the red-green-blue loops being more tightly 
nested within one another. 
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